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ABSTRACT 



Aims. To understand the origin of the Lycf line profile of the prototypical Lyman break galaxy (LBG) MS 1512-cB58 and other 
similar objects. To attempt a consistent fit of Lya and other UV properties of this galaxy. To understand the main physical parameter(s) 
responsible for the large variation of Lya line strengths and profiles observed in LBGs. 

Methods. Using our 3D Lyo- radiation transfer code (Verhamme et al. 2006), we compute the radiation transfer of Lyo- and UV 
continuum photons including dust. Observational constraints on the neutral gas (column density, kinematics, etc.) are taken from 
other analysis of this object. 

Results. The observed Lyo- profile of MS 1512-cB58 is reproduced for the first time taking radiation transfer and all observational 
constraints into account. The observed absorption profile is found to result naturally from the observed amount of dust and the 
relatively high HI column density A'h. Radiation transfer effects and suppresion by dust transform a strong intrinsic Lyor emission with 
EW(LyQr) ^ 60 A into the observed faint superposed Lyo- emission peak. We propose that the vast majority of LBGs have intrinsically 
EW(Lyci') ~ 60-80 A or larger, and that the main physical parameter responsible for the observed variety of Ly cf strengths and profiles 
in LBGs is Nn and the accompanying variation of the dust content. Observed EW(LyQ-) distributions, Lyce luminosity functions, 
and related quantities must therefore be corrected for radiation transfer and dust effects. The implications from our scenario on the 
duty-cycle of Lyo- emitters are also discussed. 

Key words. Galaxies: starburst - Galaxies: ISM - Galaxies: high-redshift - Ultraviolet: galaxies - Radiative transfer - Line: profiles 



1. Introduction 

Discovered serendipitously by lYee et al] (Il996h . the Lyman 
Break Galaxy (LBG) MS 1512-cB58 (hereafter cB58) lo- 
cated at z ~ 2.73 is strongly gravitationally magnified by the 
foreground cluster MS 1512H-36. It is by far the best stud- 
ied LBG, which has been observed at many wavelengths (see 



Baker etal. 2001, 2004; Bechtold et al. 1997; Ellingson et al. 


1996; Fraveretal. 1997; 


Nakanishi et al. 


1997; Pettinietal. 


2002, 2000; SavagHo et al. 2002; Sawicki 


2001; Teplitz et al.l 



2004.I200C 



200C ) optic al dTephtz et al.ll2004l2000h . and radio (Baker et all 
1200 IL 12004 observations have allowed to draw a fairly global 



picture of this star-forming galaxy and its interstellar medium 
(ISM): it is young (~ 50-300 Myr), massive (~ 10'° Mq), with 
a sub-solar metallicity (Z ~ 2/5 Zq), a star formation rate (SFR 
~ 24 Mq yr ') typical of LBGs, and it is surrounded by outflow- 
ing material at Vexp ~ 250 km s"' . 

Thanks to its apparent brightness the rest-UV spectra of this 
galaxy are of comparable quality as the best studied local star- 
bursts. This has in particular made it possible to detect numerous 
stellar and interstellar (IS) lines, to study their kinematics, and 
to use them for detailed UV line fits to co nstrain its stellar con- 
tent, age, star formation history etc. (e. g. Ide Mello et a n i2000t 
lElIingson et"anil996l: IPettini et alj|2000l) . Apart from its bright- 
ness, magnified b y a factor ~ 30 thanks to gravitational lensing 
jSeitz et al.ll 19981) . the properties of cB58 are empirically found 
to be typical of LBGs ( Shapley et al. 2003 ). More precisely, in 
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the sample of lShapley et al.l (l2003h cB58 is found in the quartile 
(~ 25 %) of LBGs with the strongest Lya absorption. As such a 
representative and given the quality and amout of observational 
data available, cB58 is of particular interest. 

Our objectives are to examine if the observed Lya line pro- 
file and strength of cB58 can be understood, and if one is able 
to obtain consistent diagnostics from the Lyor line and from the 
other UV and broad band features used in earlier papers. Such an 
attempt is c all ied out here fo r the first time with a Lyo' radiation 
transfer code d Verhamme et al. 2006) . Other LBGs with varying 
strengths of Lya emissi on will be analysed in a subsequent paper 
(IVerhamme et al.ll2007h . In particular we wish to elucidate what 
causes the broad Lyor absorption and very weak superposed Lya 
emissio n in cB58 and in oth er LBGs showing similar Lya spec- 
tra (cf. I Shaplev et al.ll2003l) . Similarly we want to understand 
what physical parameter(s) distinguish the different subtypes of 
LBGs showing large variations especially in Lya. Finally, fitting 
observed starburst spectra with detailed Lya radiation transfer 
modeling should in general provide insight for the use of Lya as 
a diagnostic in distant galaxies. 

The remainder of the paper is structured as follows. The main 
observational properties of cB58 are summarised in Sect.|2] The 
principles of our simulations and fit method are described in 
Sect. |3] In Sect. |4]we discuss the implications from our results 
on the understanding of Lya emission in LBGs, and on the prop- 
erties of these galaxies. Our main conclusion are summarised in 
Sect.|6] 
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2. Observational constraints from cB58 

The main properties of cB58 of relevance for the present paper 
are summarised in what follows. 



a young age is the undera bundance of N and the Fe-peak ele- 
ments in the ISM of cB58 ( Pettini et alll2002l) . giving the upper 
limit age ^ 300 Myr, the timescale for the release of N from 
intermediate-mass stars. 



2.1. Redshift and ISM outflow properties 

Direc t observations of photospheric lines bv lPettini et alj (l2002i 
20001) assign the systemic redshift at zuv = 2.7276+°°°°^. 
Teplitz et al] ( 12004) determine zhii = 2.729 + 0.001 from rest- 
frame optical emission lines(HQ;, [N ii], [O iii] and H/3) of the 
counterarc "A2", in agreement with their earlier determination 
from the main irnage of cB58. The CO emission line determina- 
tion by H^eLSLaLldZOOj) is zco = 2.7265 

The interstellar absorption lines and Lya absorption are 
blueshifted by ~ 250 km s"' with respect to the galaxy redshift, 
determined by its photospheric lines. The small Lya emission 
peak on top o f the absorpti on (see e.g. Fig. |4]i is redshifted by 
~ 300 km s ' dPettiniet al.l l2002). Both are evidence for large- 
scale outflowing medium around the star-forming galaxy. 

Interestingly, and i n contrast to th e prope rties seen in the av- 
erage LBG spectra of IShaplev et al.1 ( 2003h and spectr a of in- 
dividual z ~ 3 LBGs observed by iTapken etaH (l2O07h . Lya IS 
redshifted in cB58 by a similar amount as the interstellar lines 
are blueshifted, instead of being redshifted by — 1-2 x Vexp with 
respect to the st ellar hnes (or by ~ 3 x Vexp from the IS lines). 
As discussed by IVerhamme et al.l (1200 6) the latter behaviour is 
naturally explained by a spherically symmetric expanding shell 
with constant velocity, where the emerging Lya photons mostly 
originate from scattering on the backside of the shell where they 
gain redshift of ~ 2 x Vexp. The smaller redshift of Lya indicates 
already empirically a deviation from this simple symmetry. This 
will be substantiated below. How typical/frequent such velocity 
shifts are cannot be asserted from the available data, represent- 
ing mostly stacked spectra. This will require many high quality 
spectra of individual LBGs. However, only the Lya shape, not 
line shifts, is the criterion used to classify cB58 in the quartile 
with Lya absorption. 

The low ionisation ISM absorp tion lines yield con sistently a 
velocity dispersion b - 10 km s~' dPettini et al.ll2b02h . The neu- 
tral hydrogen column density of the outflow is determined by 
Voigt fitting of the Lya absorption wings; A^h = (7.0 + 1.5) x 10^ " 



cm ^ and A^h 



1.4 ^ jq2o 2 jg obtained by iPettini et al.l 



(I2OO2I) and ISavaglio et al.l (|2002|) respectively. Absorption pro- 
file line fitting shows that at least ~ 2/3 of the low ionisation 
materia l is found at V = - 265 km s"', with a b parameter of 70 
km S-' (iPettini et al.ll2002l) . 

The outflow of cB58 is located well in front of the stars 
and covers them almost completely, since it absorbs almost all 
the UV light from the background stars. Indeed. ISavagho et al.l 
(I2OO2I) find only a small residual mean flux above the zero level 
in the core of the Lya absorption feature, w hereas it is black 
for lPettini eFaLl (l2002l) . lHeckman et all (1200 Ih estimate an ai-eal 
covering factor for optically thick gas of 98% from the residual 
intensity at the core of the C 11 /il335 line. 

2.2. Age and star-formation history of cB58 

UV spectral fitting, including detailed stellar line profile fits, 
consistently indicate ages of ^ 20-100 Myr and favour con- 
stant star formation over this timescale (! de Mello et alj l2000t 
[Ellingson et al. 1996; Pettini et al. 2000). Such an age is also 
consistent with the observed stellar mass and star formation rate 
dBaker et al.ll2004l) . Another independent argument in favour of 



2.3. Extinction/dust content 

The precise amoun t of extinction in cB58 remains unclear. While 
IPettini et alJ (l2000l) estimated E JB - V ) = 0.29 from their UV 
spectral fitting, and lTeplitz e t al. (2000') calculated Eg{B - V) = 
0.27 from the H^/Hyj ratio, [Baker et al. (2001) used four meth- 
ods including a correction of the value derived from the Balmer 
decrement obtaining values from Es{B - V) = 0.042 to 0.236. 
Here £5 and Eg refer to the color excess of the stellar light and 
the nebular gas respectively, which are related by Es{B - V) = 
O.A4Eg(B - V) if the Calzetti law (Calzetti et al. 2000) applies. 
From this we conclude that Eg{B -V) ~ 0. 10-0.54, and we adopt 
a mean value of Eg{B -V)~ 0.3 for simplicity. 

2.4. Star formation rate, mass, and metallicity of cB58 

The star formation rate esti mated from the UV luminosity L1500 
is SFRuv ~ 40 Mq yr ' ( Pettini et al1l2000l) after correcting 
for reddening and gravitational magnification (assuming a mag- 
nification factor jj = 30, cf. below). Adopting a lower extinc- 
tion, based on the Balmer decrement. iBaker et al. I (l2004h derive 
SFRhjh^ = 23.9 + 0.7Moyr-i from Ha. 

The velocity dispersion measured from nebular and CO 
emission lines are similar: c nn - 81 km s ~' (Teplitz et al. 200Ci), 
and o-Qo = 74+18 km s ' (iBaker et al.ll20041) . The dynamical 
mass deduced fro m these values, giv en the physical size of 2.4 
kpc of the galaxy (ISeitz et al.1119981) . is Mdy„ ~ 1. x 10'° Mq. 
The mass of the gas deduced from the CO column density is 
Mgas = (6.6±^j) X 10'^ M0(Baker et al. 2003). 

The metallicity derived from the interstellar absorption lines 
and from the ionised region are co mpatible: Zis ~ 2 /5 Zq 
(IPettini et al.ll20 02) andZmi ~ 1/3 Zq dTephtz et al.ll2000l) . 



2.5. Gravitational lensing 

Lensing reconstruction shows that the primary image of cB58, 
an arc whose spectrum will be used here, is magnified by 3.35^ 
mag overall, i.e. a magnification factor fi ~ 22-AQ (iSeitz et al.l 
Il998l) . Approximately only a fraction 0.5 -0.66 of the source 
is imaged i nto th e cB58 arc according to iBaker et al.l d2004l) ; 
ISeitz etal.1 dl998l) . However, the fact that cB58 foUows all the 
same scaling relations (among Lya equiv alent width, int e rstella r 
absorption line strengths, extinction etc. IShaplev et al.l d2003l) ) 
obeyed by unlensed LBGs, shows that the spectrum is not an 
unrepresentative slice of the galaxy. 



3. Fitting the Lya spectrum of cB58 

3.1. Observations 

We use the normalised rest-frame UV spectrum from 
IPettini et al. (2002), obtained with the Echelle Spectrograph and 
Imager (ESI) on the Keck telescope II at a resolution of 58 km 
s and kindly made available to us by Max Pettini. We have 
also examined the U VES spectrum taken with the VLT from 
'Savaglio et al.l d2002h . Except for slight differences in the nor- 
malisation the spectra show a very good agreement (cf. Fig.|2]i. 
The former spectrum is used subsequently for our line profile 
fitting. 
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Table 1. Input parameters of the "standard" cB58 model for the 
radiation transfer code 



Vexp- V'front 
^back 

b 

E,(B-V) 

FWHMCLya) 

EW(LyQ-) 



255 km s^' 
free 

7.0 X 1020 

70 km 

0.3 

80 km 

free 



Vexp, and 



3.2. Lya plus UV continuum modeling 

To fit the observations we u se our 3D Monte Carlo (MC) ra- 
diation transfer code MCLya jVerhamme et alj|2006 ). The code 
solves the transfer of Lya line and adjacent continuum photons 
including the detailed processes of Lya line scattering, dust scat- 
tering and dust absorption. 

For simplicity, and given empirical evidence in favour of a 
fairly simple geome try in z ~ 3 Lyma n Break Galaxies (hereafter 
LEG) discussed by IVerhamme et al.l (l2006l) . we first adopted a 
simple "super-bubble" model to attempt fitting the observed Lya 
line profile. The assumed "standard" geometry is that of an ex- 
panding, spherical, homogeneous, and isothermal shell of neu- 
tral hydrogen surrounding a central starburst emitting a UV con- 
tinuum plus Lya recombination line radiation from its associated 
H II region. The homogeneity and a covering factor near unity 
are also supported by the observations (cf. Sect.|2]i. Taking into 
account deviations from a constant expansion velocity is neces- 
sary to reproduce the observed Lyes' profile of cB58, as already 
indicated from the empirical evidence discussed earlier. 



3.2.1 . Radiation transfer modeling 

The following geometries have been adopted: 1) a spherically 
symmetric expanding shell with velocity Vexp, 2j a foreground 
slab moving towards the observer with velocity Vfront 
3) two parallel slabs with different velocities Vback and Vfront 
Vexp respectively to mimic a shell like geometry allowing for 
velocity variations between the front and the back of the shell. 
The value of Vexp is taken from the observations. Vback is a free 
parameter. 

The other model parameters are: the velocity dispersion 
given by the Doppler parameter b, the H i column density A^h, 
and the dust absorption optical depth Tg at Lya wavelength. 
We assu me that dust and H i are uniformly mixed. As dis- 
cussed in I Verhamme et all (l2006l) . Ta relates to the usual extinc- 
tion E{B - V) ~ (0.06 .. .0.11) Ta, where the numerica l coeffi - 
cient covers a ttenuation/extinction laws of lCalzetti et al.l (l2000h . 
ISeatonI (1 19791) and similar. Here we assume E(B - V) = O.Ito. 
The values for these parameters, summarised in Table [1] are 
taken from the observations of cB58 (Sect.|2]l. Note that we sup- 
pose that the color excess measured for the nebular gas is repre- 
sentative for the dust optical depth seen by the Lya photons. 

The spectrum of the central point source is described by the 
stellar continuum around Lya, either a flat continuum or syn- 
thetic starburst spectra described below, and the nebular recom- 
bination line, described by its equivalent width ^^(Lya) and 
full width at half maximum FWHM(Lya). Although in principle 
the equivalent width of the intrinsic Lya emission is determined 
by the stellar population, i.e. can be predicted by the synthetic 
spectra, ^^(Lya) is kept as a free input parameter for the radia- 
tion transfer modeling. The resulting fit values will subsequently 
be compared to values expected from synthesis models to exam- 



ine if consistent solutions can be obtained. The value for FWHM 
is taken from the observed Ha line (cthh)- 

For the spherical shell (1), the emergent spectrum is com- 
puted by integration over the entire volume. This assumes 
that the shell is unresolved and completely integrated within 
the spectroscopic aperture, which corresponds to ~ 2 kpc 
jSeitz et al.lll998l) . For the slabs (2), the emergent spectra are 
computed taking the photons emitted within an opening angle 
of typically 60 degrees into account. The results depend little on 
the exact opening angle. 

All models have been computed using the Monte Carlo (MC) 
radiation transfer code MCLya developed by Verham me et al.l 
(2006). The code treats the scattering of line and continuum pho- 
tons in the Lya line and the adjacent continuum, as well as scat- 
tering and absorption by dust. In practice one high statistics MC 
simulation is run for each set of parameters listed above assum- 
ing a constant (flat) input spectrum. By keeping trace of the re- 
lation between the input frequency bin of the photons and their 
emergent frequency, simulations for arbitrary input spectra can 
be generated a posteriori in an interactive manner This avoids 
unnecessary transfer computations without resorting to any sim- 
plifying assumption. For comparison with the observations, the 
theoretical spectrum is finally convolved with a gaussian to the 
instrumental resolution. 

Simulations for multi-slab geometries (3) have been com- 
puted by combining the results from a single slab MC computa- 
tions in a numerical fashion, modifying appropriately the relative 
velocity shifts of the slabs and opening angles of the emerging 
radiation. This allows fast, interacti ve fitting, without res orting 
e.g. to approximate analytic fits (cf. lHansen and Ohll2006h . 

3.2.2. Syntlietic starburst spectra in tlie Lya region 

When spatially integrated galaxy spectra are considered, it is a 
priori necessary to include both stellar and nebular lines in fits of 
observed spectra. As the region around Lya is generally not con- 
sidered and not included or discussed in most evolutionary syn- 
thesis codes, we here briefly examine the behaviour of this spec- 
tr al rang e. To do so we use the evolutionary synthesis models 
of lSchaer er (2003 ) (hereafter S0 3) complemented by the recent 
models of Delgad o et al.l (l2005h including high spectral resolu- 
tion stellar atmospheres, kindly made available to us by Miguel 
Cervino. Only the main ingredients of these model sets shall be 
summarised here. 

The S03 models basically include simplified non-LTE model 
atmospheres for massive stars (> 20 Mp ), and LTE line blan- 
keted Kurucz atmospheres otherwise. The lDelgado et alj (l2005h 
models include in particular theoretical high spectral resolution 
libraries, and all massive stars are tr eated with the plane para llel 
non-LTE blanketed TLUSTY code dLanz and Hubenvll2003l) . In 
both cases the same stellar tracks and metallicities, summarised 
in S03, have been used. With respect to the Lya region these 
models represent the most up-to-date computations. The only re- 
striction is that these computations do not treat the stellar winds 
of the hottest stars, which leads to some modifications of the in- 
trinsic Lya line profile of these stars. However, unpublished test 
computations (Martins 2000) show that these effects are small 
or negligible compared to the underlying Lya absorption from 
the bulk of B and A stars and comp ared to the expec ted nebular 
emission. This is also confirmed bv lLeithereii (l2007h . 

In Fig. [1] we show the spectra predicted by these two syn- 
thesis over the spectral range centered on Lya. The spectra, 
computed for instantaneous bursts with ages 0.1, 10, and 50 
Myr at solar metallicity and for a Salpeter IMF with an up- 
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1200 1250 
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Fig. 1. Predicted stellar spectra around Lya from the synthesis 
models of S03 (blue and r ed low-resolution c urves), and the 
high-resolution models of (iDelgado et al. 1 120051) in black and 
green. Computations for solar metallicity bursts of ages 0.1, 10, 
and 50 Myr are shown from top to bottom. The green line shows 
the equilibrium spectrum reached for constant star formation, 
as the red dashed line for the S03 model. All spectra are in F,i 
units, normalised arbitrarily at 1280 A. The main photospheric 
line identifications are indicated at the top (cf. text). Note the 
excellent overall agreement of the two independent models, al- 
lowing a reasonable description of the broad absorption around 
Lya. 



per mass cut-off of 100 Mq, are normalised arbitrarily at 1280 
A. The high-resolution spectra allow the identification of the 
main (photospheric) spectral lines, which are the C iii multiple! 
centered at 1175.6 A, Si ii multiplets at AA 1193.3,1194.5 and 
1260.4,1264.7, the Si iii multiplet AA 1206.5,1207.5, and C iii 
^1247.4. The position of the N v multiplet AA 1238.8,1242.8 
formed in the winds of early type stars (not included in the 
present models) is also indicated. 

The temporal behaviour of the Lya profile in a burst can be 
sketched as follows: for ages ^10 Myr the Lya absorption deep- 
ens but remains relatively narrow. At older ages, when stars with 
a lower ionisation degree in their atmospheres dominate, the Lya 
becomes considerably stronger and broader. At > 50 Myr, not 
shown here, the red wing remains broad, whereas the flux short- 
ward of Lya disappears very rapidly. For constant (ongoing) star 
formation the spectrum reaches quite shortly (over ^10 Myr) an 
equilibrium with a Lyes' line profile quite closely resembling that 
of a ~ 10 Myr burst, as also shown in Fig.[T] 

The shape of the red wing of Lya, up to ~ 1240 A, pro- 
vides a useful constraint on the age and star formation history of 
the starburst. For bursts with ages >> 10 Myr the width of the 
stellar Lya absorption become so large, that a measurement of 
the shape of its red wing can already be distinctive, even with- 
out considering interstellar Lya radiation transfer. For example, 
for the case of cB58 shown in Fig.|2]the observed Lya wing is 
clearly less broad than the burst model at 50 Myr. Younger bursts 




-4000 -2000 2000 4000 

velocity (km/s) 



Fig. 2. Obs erved spectra of cB5 8 from iPettini et alj (|2002|) 
(black) and ISavaglio et al.l ( |2002|) (red) compared to the high- 
resolution synthesis models of jPelgado et al.l l2005l) for solar 
metallicity bursts of ages 0. 1, 10, and 50 Myr (in magenta, green 
and blue from top to bottom). Apart from small differences in 
the normalisation the two observed spectra are very consistent. 
The comparison with the models illustrates how the shape of the 
red wing of Lya can be clearly used to exclude bursts of ages 
larger than 10-50 Myr. 

or a more constant star formation history are therefore required. 
Of course, other UV spectral lines such as the well known Si iv 
and C IV lines, are se nsitive to age and star formation history (cf. 
iLeitherer etaT]|l995b . 

As other parts of the UV spectrum (see e.g. ILeitherer et al.l 

11999 ), the detailed integrated spectrum in the Lya region de- 
pends also somewhat on metallicity (cf. Fig. [3]). The main reason 
is the systematic shift of the average stellar effective temperature 
with metallicity for a given age. In consequence similar spec- 
tra of bursts can be obtained at different metallicities, provided 
some age shift is allowed for. Despite the detailed metallicity 
differences, which can be accounted for by using the appropriate 
high resolution spectra, the overall behaviour of the shape of the 
Lyof profile remains quite independently of metallicity, as could 
be expected. 

It is also reassuring to note the good overall agreement of 
the two independent synthesis models, allowing a reasonable 
description of the broad absorption around Lya even with low- 
resolution spectra. In any case predictions from both synthesis 
models, as well as flat spectra, will be used below as input for 
the cB58 fits with the Lyor radiation transfer code. 

The relative strengths of the stellar Lya absorption and ex- 
pected intrinsic nebular emission is illustrated in Fig. [3] where 
the predicted Lya equivalent widths of these components are 
plotted as a function of time for the case of instantaneous bursts 
and for constant star formation, as well as for three different 
metallicities. The behaviour of these quantities is as expected 
from earlier computations at solar metallicity ( Chariot and Fal| 
1 1 993I : fV^Us-Gabaudil 1 993b . Decreasing metallicity increases the 
nebular emission and introduces a small age shift in Lyo' absorp- 
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Fig. 3. Predicted Lya equivalent widths of the nebular emission 
(£'W(LyQ') > 0) and stellar absorption (^^(Lya) < 0) for in- 
stantaneous bursts (solid and long-dashed lines) and models with 
constant star formation (short dashed, dotted, and dash-dotted). 
The long-dashed lines show the total EW (nebular + stellar) for 
constant star formation, the dotted lines the stellar absorption 
component. Three different metallicities, Z - 0.02 (solar, black), 
0.004 (red), and 0.0004 (blue) are shown. 

tion. For constant star formation over ~ 20-100 Myr or more, 
one expects in particular ^^(Lya) ~ 60-100 A. 

4. Results 

We now discuss the simulated spectra assuming different geome- 
tries and their dependence on the input parameters. 

4.1. Models for a constant velocity shell 

The predicted Lya line profile for a symmetric spherical shell 
with constant velocity, computed with the input parameters sum- 
marised in Table [T]and with EW(LyQ')=0, 10, and 40 A is com- 
pared to the observed profile in Fig.|4] The input spectrum, here 

a flat continuum plus line, is also shown for illu stration. 

As expected from our earlier simulations (IVerhamme et al.l 

|2006|) . the global shape of the theoretic spectrum is similar to 
the observed one, consisting of a broad damped absorption plus a 
redshifted asymmetric emission line component, whose strength 
depends on the strength (EW) of the input line. The net absorp- 
tion is due to dust absorption, since the case considered here is 
the integrated spectrum of a shell, where all scattered photons 
are observed. Scattering of photons out of the line of sight is 
therefore not possible, hence any net absorption must be due to 
dust. The superposed emission peak results from the radiation 
transfer of the photons from the original (input) emission line, 
which are scattered and absorbed, and are able to escape only 
thanks to their redshift acquired on the back of the shell. 

Two main problems are immediately clear from Fig. |4] First 
the predicted width of the absorption is unsufficient, and sec- 




-4000 -2000 2000 4000 



velocity (km/s) 

Fig. 4. Predicted Lya line profile for a symmetric spherical shell 
with constant velocity, the input parameters summarised in Table 
[T] and EW(Lya')=0, 10, 40 A (green, blue, magenta respec- 
tively). The input spect rum is the red dotte d line. The observed 
spectrum of cB58 from lPettini et al.l ( |2002|) is shown in black. 



1 1 1 






1 111 ii 



-4000 -2000 2000 4000 



velocity (km/s) 

Fig. 5. Observed (black) and predicted Lya line profile for the 
single foreground slab model (blue). The input spectrum, shown 
by the red dotted line has EW(Lya)=40 A. 

ond the position of the predicted Lya emission peak is found 
at too high positive velocities. The first problem can be solved 
by increasing the column density of the shell. As already men- 
tioned above, and first pointed out by Verhamme et al. (200<^, 
an integrated shell geometry necessarily requires a larger A^h to 
reproduce a simple Voigt profile of the same width as a slab of 
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the same column density. The second discrepancy with the sim- 
ple shell model could also be anticipated (cf. Sect. |2]i, since the 
redshifted peak of Lya is expected at 2 x Vexp, whereas the ob- 
servations show the peak around ~ 300 km s"' . As the velocity 
shift of the Lya peak is a fairly generic feature of symmetric 
expanding shell models, it is clear that this discrepancy cannot 
be resolved by changes of the input parameters, such as A^h, b, 
E(B-V), etc. This has been confirmed by numerous model cal- 
culations. 

Note also that this model allows only for a fairly weak in- 
trinsic Lya emission with EW(Lya') < 10 A, much lower than 
the value expected from synthesis models for the star formation 
history derived for cB58 (cf. below). 

4.2. Single slab model 

Given the difficulties with a simple shell model and the fact that 
the direct observational constraints on the outflow only relate 
to cold material in front of the UV source, one could ask the 
question whether a single moving foreground slab (without any 
backscattering source) could explain the observations. In other 
words this boils down to the question if the observed velocity 
shift of ~ 255 km s"' between stellar and interstellar absorption 
lines could be enough to allow for a partial transmission of the 
intrinsic Lya line, such that it appears at velocities ^ 250 km 
s"'? The simple answer obtained from many simulations is no. 
More precisely, with the observed column density, b parameter, 
E(B-V), and FWHM(LyQ') the transmission of Lya line photons 
emitted at line center (V - 0) through a slab is negligible, and 
hence no line is formed on top of the broad absorption line pro- 
file, independently of the intrinsic Lya equivalent width. This is 
shown in Fig.|5] 

For enough line photons to escape trough the slab the intrin- 
sic emission line must be redshifted by > 200 km s ', which 
is larger than the possible redshift of ~ 100 km s ' found by 
iTeplitz et al.l ('2000) from 10 restfr ame optical nebular emission 
linear 

In any case, such a solution would require a large intrinsic 
Lyo- equivalent width (EW(LyQ') ^ 40 A) in agreement with the 
results found below. 

Hence the presence and position of the Lya peak together 
with a relatively large HI column density testified by the damped 
absorption excludes most likely a pure foreground slab. To allow 
for the Lyo' transmission, a receding background screen must be 
present so that the Lya photons can gain an additional redshift 
by bouncing off the back, in a similar fashion as in the spherical 
shell. Such a situation will be simulated now. 

4.3. Models for two moving slabs 

To attempt reproducing the observed Lya emission component 
we now assume a geometry with two parallel slabs, one in front 
and one in the back of the UV and Lya source. While the velocity 
and the other parameters of the approaching foreground slab are 
fixed by the observations (cf. Table [TJ, the velocity Vback of the 
receding slab is kept as a free parameters. 

A good fit to the observed spectrum of cB58 is found for 
Vback ~ 140 km s"', as shown in Fig.|6] For this choice of Vback 
the velocity shift of the receding slab is suflicient to allow a 
non-zero transmission of the reflected Lya photons through the 



' As pointed out by iPettini et alJ ( l2002h this possible small reshift 
could well be due to uncertainties in the absolute wavelength calibra- 
tion. 




-4000 -2000 2000 4000 



velocity (km/s) 

Fig. 6. Observed (black) and predicted Lya line profile for the 
two slab model described in Sect. |4.3l (blue). The input spectrum 
is shown by the red dotted line. 




-4000 -2000 2000 4000 

velocity (km/s) 



Fig. 7. Same as Fig. |6]but using high resolution spectra from 
evolutionary synthesis models for the description of the stellar 
continuum. The input spectrum, arbitrarily normalised, is shown 
by the red dotted line. 



foreground slab, and low enough to produce the peak at the ob- 
served velocity of ~ 300 km s"'. The emergent spectrum is basi- 
cally independent on the properties (e.g. A^h) of the background 
slab, which merely acts as a moving mirror This is clear since 
the spectrum reflected from a slab depends only on the ther- 
mal width b, and only weakly so (oc b^^^^), as e.g. shown by 
[Hansen and OhI (l2006l) . 
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Interestingly the value of the second free parameter, the in- 
trinsic Lya equivalent width EW(LyQ'), is required to be EW ^ 
80 A to reproduce the observed strength of the Lya emission 
peak. In our two slab model the solid angle under which the 
receding "mirror" is seen by the foreground slab governs the re- 
quired EW, since this controls the relative amount of photons 
received from direct and reflected radiation at the back of the 
foreground slab. The lower limit on EW(LyQ') given above is 
obtained by assuming a solid angle of In. The observations are 
compatible with arbitrarily larger values of EW(Lya). 

Instead of the flat continuum we have also examined models 
using the more realistic synthetic spectra including stellar Lya 
absorption presented above (Sect. 13.2.21 ) plus nebular Lya emis- 
sion as input for the radiation transfer calculations. The resulting 
Lya profile, shown in Fig. |2l is of similar quality for a slightly 
lower input Lya equivalent width of ^ 60 A. 

Fits using different combinations of the column density and 
extinction could be examined, given the uncertainty on the dust 
optical depth discussed above (cf. Sect.|2]). Using e.g. lower val- 
ues of Eg(B - V) and our "standard" value for A^h leads to a 
narrower Lya absorption profile, which would require a higher 
column density. Quantifying the degeneracies around the solu- 
tion found here is beyond the scope of this paper In any case, 
it is quite clear that the main conclusion, a high column density 
and the need for a large intrinsic Lya equivalent width, must be 
quite robust with respect to these degeneracies. This must be the 
case since a high column density is required to create the broad 
Lya absorption (independently of the exact geometry), and in 
this case the escape fraction of Lya photons "injected" close to 
line center is low quite independently of the exact dust content. 

In summary, we conclude that with the simple geometry con- 
sidered here we are able to reproduce for the first time the ob- 
served Lya profile of cB58 with a radiation transfer model us- 
ing all observed constraints on the outflow plus just two free 
parameters, Vback and EW(Lya). One of the interesting conclu- 
sions from this analysis is to show that the observed absorption 
dominated Lya line profile of cB58 is compatible with an intrin- 
sically very different spectrum with a strong (EW(Lya) <: 60 A) 
Lya emission line. Radiation transfer effects and the presence of 
dust transform the intrinsic starburst spectrum to the observed 
one. 



5. Discussion: a unifying scenario for Lya emission 
in LBGs 

As we have shown above the observed Lya spectrum of cB58 
is compatible with an intrinsically strong Lya emission (with 
EW(Lya) > 60 A). Hence this is compatible with the equiva- 
lent width expected for a starburst with a constant star forma- 
tion rate over > 50-100 Myr, the timescale over which EW(Lya) 
reaches a constant equilibrium value of ~ 60-80 A for a "stan- 
dard" IM F and m etallicities comparable to that of cB58 (see 
Fi g. [3] and ISchaerer. 2003 ) . Therefore the observed Lya strength 
and profile are compatible with the star-formation histories de- 
rived from detailed modeling of the rest-UV properties, who 
consistently found a constant star formation rate and ages up to 
ages of- 20-100 Myr (c f. [de Mello et al.ll20()0l: lElTingson etld] 
Il996t iPettini et al.l |2000 |). Constant star formation over such 
timescales, instead of an aged burst which would be responsible 
for a low EW(Lya), is also physically more plausible for objects 
like cB58 having large star formation rates, large physical scales, 
etc. Taking into account Lya radiation transfer and dust allows 
us to reconcile these diagnostics. 



Both detailed analysis of the rest-UV spectra of some LBGs 
and broad band SED fitting of many LBGs yields typically 
stellar ages of several hundred Myr a nd fairly constant or 
slowly decreasing star formation rates (cf . iPapovich et all2001l : 
iPentericci et al.l2007l : IShaplev et al.l2()0lh . These results and our 
analysis of cB58 indicate fairly lo ng star formation timescal es, 
in contrast e.g. to the suggestion of iFerrara and Ricottil (|2006|) . 

cB58 is part of the quartile of z ~ 3 LBGs showing a 
Lyg dominated by absorption with faint or little Lya emission 
dShaplev et al.l 12003 ). Therefore we can expect that these ob- 
jects behave in a similar way as cB58 and hence have much 
higher intrinsic Lya equivalent widths and higher intrinsic Lya 
luminosities. In other words observed distribution functions of 
EW(Lya) and the Lya luminosity function must be strongly 
modified by radiation transfer and dust effects, and the fraction 
of objects with low EW(Lya) (L(Lya)) must be "artificially" 
overestimated. 

Taking radiation transfer and dust into account, should there- 
fore allow to reconcile observed Lya strengths with relatively 
long (» 50 Myr) "duty cycles" of LBGs. In any case, a ge dif- 
ferenc es between ~ 200 and 400 Myr found by Penteri cci et al.l 
( l2007h for LBGs with and without Lya emission respectively, 
cannot be the direct cause of the observed Lya differences. If our 
scenario also works for Lya emitters (henceforth LAEs, i.e. Lya 
selected LBGs) and the suppresion of Lya photons is sufficient, 
it could explain the relatively low percentage of star-forming 
galaxies with no/little Lya emission and hence reduce the need 
for sho rt starburst episodes invoked e.g. bvl Malhotra and Rhoadsl 
( |2002|) . This would then also help to reconcile the lower number 
densities of LAEs com pared to LBGs despite the agreement in 
their bias jKovac et al.ll2007) . 

It may be reasonable to expect that most of the LBGs have 
quite very similar minimum intrinsic Lya equivalent widths of ~ 
60-80 A determined by constant star formation over > 50-100 
Myr Larger EWs are expected for objects dominated by younger 
($ 10^0 Myr) populations (cf. Fig. [3]). Lower EWs would then 
be expected either in post-starburst objects or due to transfer 
and dust effects. Since the outflow velocities of LBGs do not 
vary strongly between subsamples with different Lya strength^ 
we propose that the main factor leading to the diversity of Lya 
strengths in LBGs is their H i column density and concomitantly 
their dust content^. 

This suggestion is compatible with the obser vational corre- 
lations found between E(B -V) and EW(Lya) (e.g. lShaplev et al.l 
120031 : iTapken et alj 20071) and with the differences in extinc- 
tion found by IPentericci et al.l (l2007h between two subsamples 
of z ~ 4 LBGs. It also broadly agrees with Lya line profile fit- 
ting of LBGs with strong Lyg emission, discussed in the next 
paper of this series jVerhamme et alj|2007h . 

G iven observed mass- metallicity relations (e.g. lErb et al.] 

2006; Tremo nti et al.ll2004l) it may be quite natural to speculate 
that the difference in A^h and d ust content is related to the galaxy 
mass to first order. The data of lShaplev et al.l (1200 Ih show such a 
trend. In this case more massive LBGs would be more dusty, and 
hence show lower EW(Lya), although they could have identical 
Lyg emission properties intrinsically. Such a trend is e.g. found 
bv .Pentericci et al] ( l2007i) . The absence of bright LBGs (in the 



^ The expansion velocities determined by one third of the veloc- 
ity shift between the interstellar absorption lines and the Lya peak, 
VexD ~ l/3A (KisM. V(Lya)) is ~ 160-260 km s"' for the sample of 



IShapleveta i] 12003) 

^ If the dust-to-gas ratio remains constant, both A'h and the dust con 
tent increase in lockstep. 
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UV restframe) with high Lya equivalent widths (in emission) 
reported in many papers (e.g. ' Ando et al] l2004t IShaplev et al.l 
[2003; Tapken et al. 2007) would also naturally be explained by 
this fact, if the SFR correlates w ith mass, as observed for star 
forming galaxies out to z ~ 2 (see lDaddi et al.ll2007l : lElbaz et alj 
l2007tlNoeskeetal.ll2007h . 

Relevant for determining the Lya properties are however, the 
total H I column density and the relative dust to H i content. 
Assuming the latter to be constanfl, the above suggestion would 
imply larger amounts of H i in more massive LBGs, known to 
hold at least in nearby galaxies dKenni cutt 1998). However, since 
the H I responsible for shaping Lya is found in outflows it may 
only trace a fraction of the total neutral hydrogen content. Hence 
the H I column density derived from the Lya fits may not be a 
good tracer of the gas mass, and hence may not simply scale with 
galaxy mass. Correlating the outflow properties more closely to 
those of the host galaxy may require dyna mical modeling (e.g. 
iBlaizot and et aLll^OOTtTFerrara and Ricottiii2006) . 

While the "scenario" proposed here to explain the diversity 
of Lya line strengths in LBGs seems consistent with the obser- 
vations, it remains to be tested more generally. LBGs with strong 
Lya emission have also been studied with our radiation transfer 
models for this purpose and are broadly in agr eement with the 
scenario proposed here d Verhamme et al.ll2007l) . Alternative ex- 
plan ations have also been proposed. For example IShaplev et al.l 
(|2001) have from empirical evidence suggested age as the main 
difference in LBGs, where younger objects are more dusty, and 
hence show less Lya emission. However, it is not clear why older 
LBGs would contain less dust, especially since outflows, sup- 
posedly used to expel the dust, are ubiquitous in all LBGs and 
since precisely these outflows are the location where the emer- 
gent Lya spectrum is determined. Also, why would younger 
LBG s be more ma ssive than older ones as their data would im- 
ply ? iFerrara and Ricotti (2006) have proposed that LBGs host 
short-lived (30 + 5 Myr) starburst episodes, whose outflows - 
when observed at different evolutionary phases - would give rise 
to the observed correlations between IS lines and Lya. However, 
their conclusion may need to be revised since the "observed" 
wind velocity versus SFR relation (their Fig. 1) is incorrect and 
flatter than assumecfl Furthermore the ages obtained from SED 
and spectral fits of LBGs show older ages and fairly constant star 
formation histories, as alread y mentioned earlier. 

In paper III of this series (Verham me et al.ll2007h we will ex- 
amine how our models are able to explain various correlations 
observed among Ly a, IS lines, and other properties of LBGs 
(IShaplev et al.ll2003 l) and LAEs. Further work will be needed to 
test the scenario proposed here, to understand more precisely the 
relations between LBGs and LAEss, and to understand the links 
between star-formation, host galaxy, and outflow properties. 

6. Conclusions 

A 3D Lyg and UV continuum radiation transfer code 
dVerhamme et al.ll2006l) has for the first time been applied to the 



* Combining e.g. the dependence of the d ust-to-gas ratio on metal- 
licity O/H from Lisenfeld and Ferraral ( Il998l) with the mass-metallicity 
relation at z ^ 2 (Erb et al. 2006[), gives a very small dependence of the 
dust-to-gas ratio on galaxian mass. 

^ They assume v,^. = 1/2A(Vism, V(Lya)) instead of the weaker de- 
pendence Vn- ~ l/3A(VisM, V'(Lyc!')) obtained from shell models and 
supported by direct measurements of stellar, IS, and Lya velocity shifts. 
Furthe rmore their SFR values differ from those quoted in lShaplev et alj 
(I2003h . 



prototypical Lyman break galaxy MS 1512-cB58 at z - 2.7. 
Since this is one of, if not the best studied LBG, and since it 
is part of the quartile of LBGs sho wing predominantly Lya in 
absorption (cf. IShapley et al.ll2003h . a detailed Lya line profile 
analysis including radiation transfer and dust effects is of great 
interest. 

Three different geometries were explored for the material 
surrounding the central starburst; a spherically symmetric ex- 
panding shell with velocity Vexp, a foreground slab moving to- 
wards the observer, and a geometry mimicing a shell with a 
lower velocity in the back. All available observational con- 
straints were used (see Table [T]!. In particular these include: the 
velocity Vexp of the foreground material, its Doppler parameter 
b, column density A^h, and the observed extinction. The width of 
the intrinsic Lya emission was taken from the observed FWHM 
of Ha. The two free parameters were: the equivalent width of the 
intrinsic Lya emission, and the velocity Vback of the background 
slab (if applicable). 

Our radiation transfer calculations have confirmed what 
could be expected from our earlier expanding shell simulations 
and from the lower-than-usual observed velocity shift between 
the peak of "remnant" Lya emission and the blue IS absorp- 
tion lines (cf. Sect. 12.11 1; the Lya line profile of cB58 cannot 
be reproduced by an expanding shell with isotropic constant ve- 
locity Vexp. Considering, however, a lower shell velocity in the 
back allowed us to obtain excellent fits to the observed Lya line 
profile. Interestingly, the input spectrum requires a fairly strong 
Lya emission, with a lower limit of EW(Lya) ^ 60 A. Radiation 
transfer and the suppresion of photons by dust are responsible 
for transforming such an input spectrum into the observed Lya 
absorption profile with a superposed faint emission peak. In this 
way the observed Lya line strength and profile can be reconciled 
with the strong intrinsic Lya emission expected from the approx- 
imately constant star-formation history o f cB58 derived from 
earl ier detailed U V spectral analysis (cf . Ide Mello et al.] 120001: 
■Pettini et al.ll2000h . 

In fact we suggest that cB58 and most other LBGs have 
intrinsically EW(Lya) ~ 60-80 A or larger, and that the main 
physical parameter responsible for the observed variety of Lya 
strengths and profiles in LBGs is A^h and the accompanying vari- 
ation of the dust content (see Sect. [Sj. This explains not only 
the absorption-dominated object cB58, but also observed corre- 
lations between E(B-V) and EW(Lya) in LBGs, the absence of 
bright LBGs with strong Lya emission, and other correlations. 

Among the implications from our work are that observed 
EW(Lya) distributions and Lya luminosity functions must be 
corrected for radiation transfer and dust effects. Furthermore rel- 
atively short duty cycles, suggested earlier in the literature, are 
not required to explain the variations observed between different 
LBG types. Our proposed unifying scenario will be detailed fur- 
ther and subjected to additional tests in a subsequent publ ication 
on a sampe of LBGs with strong Lya emission (Verhamm e et al.l 
'2007). Providing a clearer picture of the physical links between 
the observables including Lya line strength and profile, and the 
star-formation, host galaxy, and outflow properties and evolution 
of LBGs and LAEs remains an objective for the near future. 
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